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RESUMEN 

Para el analisis no lineal de estructuras, es necesario conocer el diagrama 
momento- rotacion o el diagrama momento-curvatura, para las vigas y columnas 
de hormigon armado o de acero y para las diagonales rigidizadoras se necesita 
conocer el diagrama fuerza-desplazamiento. El ASCE 41 de 201 1 , proporciona 
tablas que facilitan el calculo de estos diagramas, las mismas que se presentan en 
este articulo, con ejemplos desarrollados, manualmente. 

Entre los ejemplos se presenta una estructura de 4 pisos de hormigon 
armado y una estructura de acero de 6 pisos; en las dos, se tiene disipadores de 
energia ADAS o TADAS, colocados sobre diagonales de acero en forma de "V" 
invertida; se presentan los diagramas momento curvatura en varies elementos de 
estas dos estructuras y manualmente se ilustra el calculo del diagrama fuerza- 
desplazamiento para las diagonales de acero, tanto para el case de compresion 
como para el case de traccion; de igual forma se obtiene manualmente el 
diagrama momento-curvatura para un disipador de energia. 

El calculo con las tablas del ASCE 41 implica evaluar algunos criterios y la 
mayor parte de las veces se debe interpolar con los valores de las tablas, de tal 
manera que es bastante laborioso hacerlo manualmente por lo que se presentan 
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varios programas del sistema de computacion CEINCI-LAB que obtienen los 
diagramas indicados en el primer parrafo. 



ABSTRACT 

For nonlinear analysis of structures, it is necessary to know the moment- 
rotation diagram or moment curvature diagram, for beams and columns of 
reinforced concrete or steel and for the stiffening diagonals, these needs to know 
the force-displacement diagram. The ASCE 41 of 2011 provides tables that 
facilitate the calculation of these diagrams, the same that those presented in this 
article, with examples developed manually. 

In the examples are presented a reinforced concrete structure of 4 floors 
and a steel structure of 6 floors, in the these structures, they have energy 
dissipators TADAS or ADAS mounted on steel bracing in the form of inverted "V"; 
the moment curvature diagrams are presented in various elements of these two 
structures and manually calculated the force-displacement diagram for steel 
diagonal, for both the case compression as in the case of traction, likewise 
manually, it is obtained the moment-curvature diagram for an energy dissipators. 

The calculation with the table of ASCE41 involves assessing certain criteria 

and most of the time must be interpolated with the values of the tables, so it is 
quite laborious to do it manually, so several programs are presented of computer 
system CEINCI-LAB. They get the diagrams indicated in the first paragraph. 

1. INTRODUCCION 

Para el analisis no lineal: estatico o dinamico, se necesita conocer las 
relaciones momento-curvatura o momento-rotacion, para el caso de flexion y flexo 
compresion; de estos diagramas se obtiene la rigidez de una seccion de acuerdo 
al nivel de cargas que gravita sobre la misma. (Aguiar et al., 2015, 1; Mora y 
Aguiar, 2015; Medina, 2012). 

En el sistema de computacion CEINCI-LAB se tienen los programas 
MCviga y MCcolumna que determinan el diagrama memento curvatura de vigas y 
columnas, respectivamente, utilizando un modelo de hormigon no confinado y un 
modelo elasto plasto para el acero. Ortiz y Freire (2013) 

Para hormigon confinado tambien se han desarrollado programas para 
encontrar el diagrama memento curvatura, utilizando el modelo tri lineal para el 
com portam lento del acero. Aguiar (2002). En resumen el diagrama memento 
curvatura depende de los modelos constitutivos del hormigon y del acero que se 
utilicen y para cualquiera de los casos con los que se trabaje se necesita contar 
con programas de ordenador. Borda (2007), Chioul et al. (2008), Rogac (2012), 
Lepe (2013). 

Ahora bien, el ASCE 41 presenta unas tablas que tienen un respaldo 
teorico y experimental, por medio de las cuales se puede encontrar en forma 
sencilla las relaciones memento rotacion y a partir de estas inferir las relaciones 
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momento curvatura con las que se halla la rigidez de acuerdo en una seccion de 
acuerdo al estado de dano que tiene la misma. 




Q/Qv 



eo A 



Figura 1 Relaciones momento curvatura y modelo del ASCE 41 



A la izquierda de la figura 1, se presenta las relaciones momento 

curvatura, para un modelo que contempla tres zonas de trabajo del elemento; la 
elastica (sin dano) que es para cuando el momento actuante M < My] siendo My 
el momento de fluencia; en esta zona la rigidez a flexion vale E Ig. Donde E es el 
modulo de elasticidad del material; 4 es el momento de inercia de la seccion 
gruesa, que se halla en base a la geometria de la seccion; asi para una viga 
rectangular de seccion constante Ig = b h^/12, donde b,h, son la base y altura de 
la seccion; de tal manera que la rigidez elastica se halla en base a la geometria de 
la seccion transversal del elemento sin utilizar el ASCE 41. De igual manera se 
halla el momento de fluencia My. 

La segunda zona es la de post fluencia, en esta se cumple My < M < 
My, y la rigidez de la seccion es E Ip] para esta zona se requiere las propuestas 
del ASCE 41 pero antes de explicar su calculo y para terminar con las zonas se 
indica que la ultima es la residual que tiene un momento M^. 

El diagrama presentado a la izquierda de la figura 1 es para el caso de 
secciones asimetricas que tienen un com portam lento diferente si la parte inferior 
esta traccionada y la superior trabaja a compresion; o al reves. 

A la derecha de la figura 1 , se presenta la nomenclatura del ASCE 41 para 

hallar las relaciones momento rotacion. En el eje horizontal se indica 9 (giro) o 
desplazamiento A; y en el eje vertical Q/Qy Para el caso de flexion esta relacion 
sera M/My y lo interesante es que el valor para el punto B, es la unidad, de tal 
manera que el momento es el de fluencia My, que esta asociado a una rotacion By 

A partir del punto My, By, se encuentra el punto U (grafico de la izquierda) 
que corresponde al punto C (grafico de la derecha) por medio de la variable a 
propuesta por ASCE 41. Se determina tambien el momento residual, M^, con la 
variable c; y la rotacion final en este segmento DE por medio de la variable b. Se 
destaca que en el hormigon el punto U esta asociado a una deformacion a la 
compresion que esta entre 0.003 y 0.004; en cambio para la zona residual esta 
deformacion se halla entre 0.004 y 0.008. 
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Las variables a,b,c, presentan el ASCE 41 de 2011 para elementos de 
acero y hormigon armado y se indican en este articulo con ejemplos y programas 
de computacion del sistema CEINCI-LAB. 

Es una gran ventaja contar con las tablas del ASCE 41 para hallar las 
relaciones momento-rotacion, en forma relativamente sencilla pero se debe tener 
claro que en el modelo del ASCE 41 (figura 1) el punto C, se encuentra a la 
derecha del punto B, es decir la seccion tiene ductilidad por curvatura pero si la 
viga no tiene ductilidad debido a que la cuantia de acero a traccion es similar o 
igual a la cuantia balanceada, no va a existir el punto C. (Marm 1979). De tal 
manera que se recomienda un poco de precaucion, tratando de no estar en ios 
limites de Ios criterios que presentan las tablas. 

2. PUNTO DE FLUENCIA PARA VIGAS Y COLUMNAS DE 
HORMIGON 



Como se acaba de ver el diagrama memento rotacion se halla a partir del 
punto de fluencia Y, utilizando Ios coeficientes a,b,c, del ASCE 41. Para encontrar 
el punto Y se decidio trabajar con la propuesta de Y. Park (1985) que cuenta con 
un respaldo teorico y experimental del ensayo de 400 elementos de vigas y 
columnas de hormigon armado. Las ecuaciones de calculo se indican a 
continuacion. 

My = 0.5 b [(1 + /?c - /?) /7o + (2 - J?) Pt+ ir]- Z/JJa^ Ptl 
d' _ 0.75 _£y _ Po 



" "d ^ " 1+ a, W 



fy , fy 



Sy 



0, = [1.05+ fe- 1.05)^]^' 



0.03J (1 - k)d 



1 11 



k= (Pt + v'tY ir^+ iVt+ Pc V't) — - (Pt + V't) 



C2 = 1 + 



0.45 



(0.84 + pt) 
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Donde b,d,d', son la base, la altura efectiva y ei recubrimiento de la 
armadura a compresion; fy.Sy.fc.i^o son: el limite de fluencia del acero, la 
deformacion de fluencia, la resistencia maxima del hormigon a compresion, y la 
deformacion asociada a la resistencia maxima del hormigon; Pq es la carga axial 
que gravita sobre el elementos, para vigas es igual a cero; Ec la deformacion del 
acero a compresion que tiene que ser menor a la deformacion ultima £„ que puede 
ser igual a 0.003 para disefio o 0.004 para analisis; A^.A'^ es el area de la 
armadura que trabaja a traccion y compresion respectivamente, para el caso de 
columnas se considera que ^4^ = A'^; My,0y. son el momento y curvatura en el 
punto de fluencia. 

Con el modelo de Park (1985) tambien se puede hallar el momento y 
curvatura para el punto ultimo (My,0„), en Aguiar (2003) se indica las ecuaciones 
de calculo para este punto. 



3. RELACIONES PARA VIGAS DE HORMIGON 

En las tres primeras columnas de la tabia 1, se indican los parametros que 
se deben evaluar, antes de encontrar los valores de a, b. c. que estan en las tres 
ultimas columnas. El primer parametro esta relacionado con las cuantias de acero; 
el segundo define si el refuerzo transversal utilizado esta confinando o no a la viga 
y el tercero esta relacionado con el cortante actuante. A continuacion se indica el 
calculo de estos parametros. 

• Cuantia a flexion 

As , A's 0.85 /c' / 6120 

P = r^ P=r3 Pb= — ? — P 



b~d ^ ~ b~d ~ fy V6120 + /y^ 

= 0.85 para/'c < 280 kg/cm} 

= 0.65 < 1.05 - /'c/1400 < 0.85 Para/'c > 280 kg/cm^ 

Las variables todavia no definidas son: p.p'.Pt,, que son la cuantia del 
acero a traccion, a compresion y la cuantia balanceada. Con estas 

cuantias de determina la relacion: y en base a su valor se ingresa a 

la tabIa 1 . 

Refuerzo Transversal confinado, si cumple: 

d 3 Ay* f-y* d 
s<- Vs>-V Vs = 

Donde A^ es el area del refuerzo transversal; fy es el limite de fluencia; d 
es la altura efectiva; s es el espaciamiento de los estribos; 6S la 
capacidad al corte de la armadura transversal; V es el cortante actuante. 
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• Cortante actuante 

V 

Si S8 trabaja el cortante 7 en T.; las dimensiones de la viga o columna: 
h„,d, en m., y la resistencia del hormigon en T/m^, se debe multiplicar el 
resultado por 1.1926 para pasar Pound/in. 



Tabia 1 Parametros y criterios de aceptacion para procedimientos no 
lineales en vigas de hormigon. 



i.i) Vigas controladas por flexion 


a 


b 


c 


P-P' 


Refuerzo 
Transversal 


V(pound) 








Pbal 


by„{in)d{in)^J f' c{psi) 


< 0.0 


Confinado 


< 3 


0.025 


0.05 


0.2 


< 0.0 


Confinado 


> 6 


0.02 


0.04 


0.2 


> 0.5 


Confinado 


< 3 


0.02 


0.04 


0.2 


> 0.5 


Confinado 


> 6 


0.015 


0.02 


0.2 


< 0.0 


No Confinado 


< 3 


0.02 


0.03 


0.2 


< 0.0 


No Confinado 


> 6 


0.01 


0.015 


0.2 


> 0.5 


No Confinado 


< 3 


0.01 


0.015 


0.2 


> 0.5 


No Confinado 


> 6 


0.005 


0.01 


0.2 


i.ii) Vigas controladas por corte 








Espaciamiento de estribos < d/2 


0.0030 


0.02 


0.2 


Espaciamiento de estribos > d/2 


0.0030 


0.01 


0.2 



• Ejempio 1 

Dada una viga cuadrada de 40 cm., de lado y doblemente armada como 
se indica en la figura 2; el recubrimiento es de 6 cm.; el refuerzo transversal esta 
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constituido por varillas de 10 mm, espaciados cada 10 cm. El hormigon tiene un 
f c — 240 kg/crn}.^ el acero es A 36. El modulo de elasticidad del hormigon es 
1500000 T/m2 , del acero 20430000 T/m^. Se pide: 

i) Determinar la rigidez axial EA, y rigidez a flexion EI en el range elastico. 

ii) Encontrar los parametros a,b,c, de la tabia 1 si el cortante 7 = 0 y el 
memento inicial M = 0, inicio de calculo. 

ill) Hallar el Memento y curvatura a nivel de fluencia. 



Ec = I'SOO.OOO T/m^ 



fy = 25300 r/m^ 
Es = 20'430.000 T/m^ 
fu = 40800 T/m^ 
Figura 2 Geometrfa de una seccion de una viga y parametros del 
acero y concrete. 



1 E01Omm ©lOcrn 



• 3020mm 
o 3016mm 



• Solucion 

1) Rigidez axial y a flexion 

A = 0.40 * 40 = 0.16m^ 

0.40 * 40^ 

/ = — = 0.00213 m* 

12 

La rigidez axial EAy a flexion EI, valen: 

EA = 240000 T. EI = 3200 T/m 

ii) Factores a, b, c, de tabIa 1 . 

As = 3016 mm = 6.03 As' = 3020 mm = 9.42 cm^ 

6.03 9.42 

p= ^ = 0.0038 p' =——-= 0.0059 p^, = 0.0053 

40*34 ^ 40*34 

p-p' 

- — ^= -0.4039 
Pb 



198 



Roberto Aguiar, David Mora y Michael Rodriguez 



A continuacion, se procede a determinar si la seccion se encuentra 
confinada, para lo cual la viga se encuentra armada con 10 10 mm @10cm. 

= 2010 mm = l.S7cm^ 

1.57*4200*34 

Vs = — = 22.43 T 

^ 10 

Al cumplir las dos condiciones presentadas anteriormente se puede 
determinar que la seccion se encuentra confinada. Por otra parte, al ser 7 = 0, 
dato del problema. Los valores a,b,c son los de la primera fila de la tabia 1. 

a = 0.025 b = 0.050 c = 0.200 

ill) Momento y curvatura de fluencia 

Al aplicar las ecuaciones de Y. Park (1985) se halla que el momento y 
curvatura de fluencia valen: 

Tdd 

My = 7.6679 Tm. 0y = 0.0081 

m 

3.1 Longitud plasties 

Cuando se supera el momento de fluencia My se tiene una longitud 
plastica en el nudo inicial que se denomina L^^, y cuando se supera el momento 
de fluencia My en el nudo final se tiene la longitud Lp^; a la izquierda de la figura 3 
se presenta la nomenclatura indicada donde los mementos actuantes se han 
denominado Mj,M^; en la grafica estos dos mementos son anti horarios. 




Figura 3 Diagrama de mementos y longitudes plasticas 



En la grafica de la derecha de la figura 3, el valor de Mj se lo ha colocado 
al lado del momento Mj, ya de que de esta forma se facilita el calculo de la longitud 
plastica. En efecto, de la relacion de triangulos semejantes, se tiene: 

+ Mj My + Mj 
L L — Lpii 
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Al despejar Lp^ se encuentra: 

^ /M; - Mv 



Para el nudo final, se procede en forma similar y se halla: 



\Mi+ Mj 

Notese en las dos ecuaciones que el numerador debe ser siempre 
positive, de tal manera que la longitud plastica sera positiva. 



3.2 Rotacion de fluencia 



Para calcular la rotacion de fluencia By es necesario recordar la forma 
como se obtiene la matriz de flexibilidad en una viga que se encuentra sobre un 
apoyo fijo en el nudo inicial y un rodillo en el nudo final. Aguiar (2014). 

Para encontrar la matriz de flexibilidad se aplica un momento unitario en el 
nudo inicial y se calculan los giros en el nudo inicial y final que son L/3EI, L/6EI, 
respectivamente, grafico superior izquierdo de figura 4; al aplicar el momento 
unitario en el nudo final se hallan los giros indicados en la grafica superior derecha 
de la figura 4. Se recuerda la convencion de signos que anti horario es positive. 




Figura 4 Elementos de la matriz de flexibilidad de un elemento de seccion 

constante 



Ahora, cuando se aplican los dos mementos unitarios a la vez, al aplicar el 
principle de superposicion lineal se tienen los giros indicados en la parte inferior de 
la figura 4. 

Si en lugar de aplicar un momento unitario, se tiene un momento My en el 
nudo inicial, y un momento My en el nudo final; los giros que se encuentran son los 
de fluencia 6y,6y- Nuevamente per el principle de superposicion, estos giros valen: 
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LMy lM!y 

3EI ~ 6 EI 

LM'y LMy 

3EI ~ 6 EI 




LM;/3EI-LM,/<5EI / mSk 
Figura 5 Rotaciones de fluencia en nudo inicial y nudo final 

La figura 5 considera que las dos secciones ingresan simultaneamente al 

rango no lineal (hipotesis conservadora) ya que la mayor parte de las veces solo 
una de las secciones ingresara ai rango no lineal y la otra no. 

Para el case particular en que My = My, la rotacion de fluencia es la 
recomendada por el ASCE 41, que resulta la siguiente para elementos de 
hormigon armado. 



ey = 



LMy 

6 EI 



Para elementos de acero el memento de fluencia My = z fy. Luego la 
ecuacion resulta. 



Las variables todavia no indicadas, son: z es el memento estatico (z = 
fydA = Ac* brazo_s + At * brazoj); fy es el esfuerzo de fluencia del acero. 

Las ecuaciones indicadas son para vigas. Ahora para columnas, el ASCE 
41 recomienda la siguiente ecuacion. 




Donde P es la carga axial que gravita en la columna; Pcl la carga critica 
de pandeo. De acuerdo al ASCE 41, la carga critica se calcula de la siguiente 
manera. 

Se determina la relacion L/r. Donde L es la longitud del elemento; r es el 
radio de giro. Se debe verificar que: 
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L 

- < 108 
r 

Si la relacion L/r es menor a 108; el esfuerzo del material F^^ es 12 (psi); 
caso contrario vale: 

1.4* 10^ 



(f) 



Para pasar de (psi) a (T/m^) se debe multiplicar por 703.069. Finalmente, 
la carga critica de pandeo se halla con la siguiente ecuacion. 

P=AF,r 

Donde A es el area de la seccion transversal. 

3.3 Factor que mide efecto de corte 

El factor que mide el efecto de corte 0 para un elemento de seccion 
constante que trabaja en el rango elastico vale (Aguiar, 2014) 

3 pEI 

Las variables no definidas son: ^ coeficiente de forma que vale 1.2 para 
secciones rectangulares; G es el modulo de corte. Ahora bien, interesa calcular el 
factor 0 para un elemento de seccion constante, que ha ingresado al rango no 
lineal y tiene una longitud plastica L^a asociada a una rigidez £■/„, en el nudo 
inicial; en el nudo final la longitud plastica es 1^,^ y la rigidez es EI^; la longitud 
central que no ha ingresado al rango no lineal es: L-Lp^- Lpb, y tiene una 
rigidez EIq. Solo para el calculo del factor 0 se determinar una rigidez ponderada 
en funcion de las longitudes, denominada 



FI = 



^pa ^^a ~ ^pa ~ ^pb) ^^0 + ^pb 



Finalmente se determina un solo factor del efecto de corte, aplicabie ai 
nudo inicial, centro de luz y nudo final, denominado 0^ 



0m = 



3 (B EI 



m 



GAL2 



• Ejempio 2 

La viga de 40/40 del Ejempio 1 , tiene una longitud de 6.0 m.; un momento 
de fluencia My = 7.6679 Tm., una curvatura de fluencia 0y = 0.0081 -. En un 

m 
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instante de cargas dinamicas, el momento actuante en sus extremos es de Mj 
Mj = 8 Tm., y el cortante actuante es V = 2.7 T. Se pide: 

i) Calcular la rotacion de fluencia dy de acuerdo al ASCE 41 ; 

ii) La longitud plastica, considerando My = My. 

iii) El diagrama Momento rotacion con la formulacion ASCE 41 

iv) El diagrama Momento curvatura de acuerdo al ASCE 41 



• Solucion 

i) Rotacion de fluencia 

6 * 7.6679 
= . = 0.0024 

^ 6* 3200 

ii) Longitud plastica 

Los momentos actuantes en los extremos son anti horarios y valen 8 Tm., 
(dato del ejempio), los momentos de fluencia valen 7.6679 Tm, como se 
indica en la figura 6. 



8.00 




0.1245 

^^ ^ 



QJ24^ 






8.00 



6.00 



Figura 6 Momentos actuantes y de fluencia para calcular longitud plastica 



/8.00 - 7.6679\ 
= Htttt^ — ^TTT^ * 6 = 0.1245 m 
V 8.00 + 8.00 / 



iii) Relacion momento rotacion 

En el Ejempio 1 , se encontro que (p - p')/Pb es menor a cero; que la viga 
esta confinada; y ahora, en base al cortante actuante de V = 2.7 T. Se 
demuestra que la relacion V/(b d JJJ) es menor a 3. 



V 



2.7 



hd^JI 0.4 * 0.34 * V240 



* 1.1926 = 1.5283 



Pound 



in 



Luego de la tabia 1 se encuentra: 
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a = 0.025 b = 0.050 c = 0.200 
Con estos valores se halla el diagrama momento rotacion de acuerdo al 
ASCE 41 

9u = 9Y + a = 0.0024 + 0.025 = 0.0274 rad 
Mu = My + 0.05 EI (Bu - Oy) 
My = 7.6679 + 0.05 * 3200 (0.0274 - 0.0024) = 11.6679 T m 
Mji = cMY = 0.20 * 7.6679 = 1.5336 T m 
eii = 6Y + b = 0.0024 + 0.05 = 0.0524 rad 
El diagrama momento rotacion se indica en la figura 7. 



12 



10 



8 



4 



2 



0 



0 0.01 0 02 0 03 0 04 0 05 0.06 

Rotacion (rad) 



Figura 7 Diagrama momento rotacion de acuerdo al ASCE 41 



iv) Diagrama momento curvatura 

rad 

0Y = 0.0081 

m 
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e, 



0.0274 



0, = 0, + ^= 0.0081 + ^= 0.2282. 



rad 
m 

rad 



9u 0.0524 
0„ = 0, + -= 0.0081 + — = 0.4290 ^ 



En la figura 8 se indica el diagrama momento curvatura. 




0 2 0 25 

Rotaclon (rad'm) 

Figura 8 Diagrama momento curvatura de acuerdo al ASCE 41 

4. RELACIONES PARA COLUMNAS DE HORMIGON 

En forma similar a las vigas, hay tres controles que se deben realizar, antes de 
ingresar a la tabia 2 para encontrar los coeficientes a,b,c con los cuales se halla el 
diagrama momento rotaclon. El primer control es calcular la relacion: 



045 



Aq fc 



Donde P es la carga axial; Ag es el area de la seccion transversal de la 
columna y fc es la resistencia maxima del hormigon a compresion. El segundo 
control tiene que ver con el confinamiento y el tercero con el cortante actuante; 
estos dos controles se realizan en forma similar a la indicada para vigas. 

Si la columna esta controlada por corte, ultima fila de la tabIa 2, no se 
puede obtener los coeficientes a,b,c. 

TabIa 2 Parametros y criterios de aceptacion para procedimientos no lineales 
en columnas de hormigon armado. 



ii.i) Columnas controladas por flexion 


a 


b 


c 
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p 


Refuerzo 
Transversal 












by^,(in)d(in)y/ f'c(psi) 


< 0.1 


Confinado 


< 3 


0.0 

2 


0.0 

3 


0. 

2 


< 0.1 


Confinado 


> 6 


0.0 

16 


0.0 

24 


0. 

2 


> 0.4 


Confinado 


< 3 


0.0 
15 


0.0 

25 


0. 

2 


> 0.4 


Confinado 


> 6 


0.0 

12 


0.0 

2 


0. 

2 


< 0.1 


No 

Confinado 


< 3 


0.0 

06 


0.0 

15 


0. 

2 


< 0.1 


No 

Confinado 


> 6 


0.0 


0.0 

iz 


0. 

z 


> 0.4 


No 

Confinado 


< 3 


0.0 

03 


0.0 

1 


0. 

2 


> 0.4 


No 

Confinado 


> 6 


0.0 
02 


0.0 
08 


0. 

2 


ii.ii) Columnas controladas por corte 








d P 

Espaciamiento de estribos< o — -7- < 0.1 

^ 2' Agf'C 


No se permitira 



• Ejempio 3 

La columna indicada en la figura 9 es cuadrada de 50 cm., de lado; armada 
con 8 varillas longitudinales de 20 mm., de diametro y el refuerzo transversal es de 

1 0 de 10 mm., espaciado cada 10 cm., incluye adicionalmente una vincha en 
cada direccion del mismo diametro. La resistencia del hormigon es f'c = 
240 kg/cm^.y el acero es A 36. El modulo de elasticidad del hormigon es 
1500000 T/m2 , del acero 20430000 T/m^ Se pide: 

i) Determinar la rigidez axial y a flexion en el rango elastico. 

ii) Si la carga axial P = 20 7 y el cortante es V = 2T. Determinar los 
coeficientes a, b, c de la tabia 2. 

ill) Detallar el calculo del momento y curvatura de fluencia. 
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1 E01Omnn @10cm 





■0 — 


J? — 


2fl S 



m 3020 mm 
« 2020 mm 
O 3020 mm 



.50 



Flgura 9 Geometria da columna da Ejampio 3. 

• Solucion 

i) Rigidez axial y a flexion 

A = 0.50 * 50 = 0.25m^ 
0.50*50^ 



/ = 



12 



= 0.0052 



La rigidez axial EAy a flexion EI, valen: 

EA = 375000 T. EI = 7812.50 T/m 

ii) Factores a, b, c, de tabia 2. 

P 20.00 

= 0.0333 



Ag fc 0.25 * 2400 



2.0 



:* 1.1926 = 0.2213 



Pound 



bd^fl 0.5 * 0.44 * V240 
Verificacion de confinam lento de la seccion 

Ay = 3010 mm = 2.35 cm^ 
2.35*4200*44 



in 



10 



= 43.43 T 



d 
3 



s < - = 10 < 14.66 Vs>-V = 43.43 > 1.5 

s 4 



Diagrama momento-rotacion y momento-curvatura 



207 



Al cumplir las condiciones de confinamiento de la seccion, se selecciona los 
valores de a, b, c, de la Tabia 2. 

a = 0.020 h = 0.030 c = 0.200 

iii) Detallar el calculo del momento y curvatura de fluencia. 

Al utilizar las ecuaciones de Y. Park (1985) se halla que el momento y 
curvatura de fluencia valen: 

Tdd 

My = 27.2074 Tm 0y = 0.010 

m 

Se calcula la Rotacion de fluencia con la formula mencionada anteriormente 

3.00* 27.2074 

= . ^^.or-^ = 0.0017 rad 
y 6*7812.50 

Con estos valores se halla el diagrama momento rotacion de acuerdo al 
ASCE 41 

Bu = 9Y + a = 0.0017 + 0.020 = 0.0217 rad 
dii = 6Y + b = 0.0017 + 0.03 = 0.0317 rad 
Mu = My + 0.05 EI (Bu - 9y) 
My = 27.2074 + 0.05 * 7812.50 * (0.0217 - 0.0017) = 35.02 T m 
Mji = cMY = 0.20 * 27.2074 = 5.4415 T m 




0 015 0 02 

Rotacion (Rad) 



0.035 



Flgura 10 Diagrama momento rotacion de acuerdo al ASCE 41 
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5. ESTRUCTURA DE HORMIGON ARMADO CON 
DISIPADORES 



Una forma eficiente de reforzar estructuras es mediante la colocacion de 
disipadores de energia, razon por la cual se han desarrollado una gran variedad 
de estos dispositivos de control pasivo. Tena (2000). 

Dentro de esta gran variedad se tienen los disipadores por fluencia del 
material y aqui se tiene otra serie de sistemas de control. (Nuzzo et al. 2015). De 
tal manera que el tema ha side ampliamente investigado por la bondad que estos 
representan para las estructuras. 

En este articulo se presentan dos tipos de disipadores de energia por 
fluencia del material y son los denominados ADAS y TADAS (Tena, 1997; 
Withaker et al., 1989; Aguiar et al., 2015,1). Los mismos que se encuentran 
colocados sobre contravientos de acero en forma de "V" invertida. Concretamente 
en este apartado se ilustrara como se halla el modelo constitutivo para los 
elementos de acero, tanto para el case de compresion como para el case de 
traccion y en el siguiente apartado se hara lo propio con los disipadores. 



A la izquierda de la figura 11, se muestra un portico piano de hormigon 
armado de 4 pisos; las dimensiones y armadura de las vigas y columnas se 
indican a la derecha de esta grafica y son las que han servido de base para la 
realizacion de los ejemplos 1 a 3. 



0.75 T/m 



0.75 T/m 






1 E01OmmeiOcm 



r 









• 3020 mm 

• 21320 mm 
031320 mm 



SEOCION A-A 
COLUMNA 



1 E01Omm91Ocm 



• 3020mm , 

o ^es^6^wn 



SECadNB-B 




SECCION D-D 
DIAGONAL 



1 E01Omm92Ocm 



• 2)320mm 
o 301<mm 



secci6nc-c 



VIGAS 
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Figura 11 Geometria de un portico piano de 4 pisos de hormigon armado, con 
disipadores de energia colocados sobre contraviento Chevron. 

En la figura 12 se indica la geometria de los disipadores de energia 
colocados, a la izquierda de esta grafica se indican los TADAS y a la derecha los 

ADAS; los disipadores del primer piso son diferentes, tienen 6 placas de acero y 
los disipadores de los pisos dos y tres son iguales, tienen 4 placas de acero. Las 
diagonales son de acero compuestas por un perfil tubular cuadrado de 100/100/10 
mm., como se indica en la parte superior derecha de la figura 1 1 . 



GEOMETRIA DE LA PLACA 
nm .0.10, 



GEOMETRIA DE LA PLACA 



111 nn 


nn iil 












Figura 12 Geometria de los disipadores TADAS a la izquierda y ADAS a la 

derecha. 



.0.10 



).01 



Lo_p L 

0.20 

i 



.0.10, 



101 



0.20 
I 



En la figura 13, se presenta a la izquierda la numeracion de los nudos y 
elementos. Notese que al disipador se modela como un elemento corto y que las 
coordenadas principales van del 1 al 7. (Grafico que esta a la derecha de la figura 
13). 
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Figura 13 Numeracion de nudos, elementos y grados de libertad. 



En la figura 14 se indica los diagramas momento-curvatura, hallados en la 

columna 1 y en la viga 13 (ver figura 13) considerando un valor de a = 0.05 
(relacion entre rigidez post fluencia con respecto a rigidez elastica) y una ductilidad 
de curvatura de fi^ = 10. 
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Elem=l 




0 0.005 0.01 0.015 0.02 0 0.005 0.01 0.015 

Curvatura(rad/m) Cufvatura(rad/m) 

Figura 14 Diagramas momento-curvatura de una columna y una viga de hormigon 

armado. 

En la figura 15 se presenta el diagrama momento curvatura hallado para 
los disipadores, en la cara superior, son iguales para el ADAS y TADAS; debido a 
que los dos en la cara superior tienen un ancho b = 10 cm., y una altura de la 
seccion transversal t = 1 cm. En base a estas dos dimensiones: 10/1 cm., se 
determine el punto de fluencia, teniendo en cuenta que es de acero el elemento y 
para el punto U. se trabajo con el mismo a = 0.05 y la misma ductilidad de 
curvatura, H0 = 10 , debido a que el objetivo es presentar ejemplos del diagrama 
momento curvatura. 



Por la geometria de los disipadores se espera que el dano se produzca en 
todo la placa de tal manera que el diagrama momento curvatura indicado en la 
figura 15 se aplica a todo el elemento disipador. 
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• Ejempio 4 

Determinar el modelo constitutivo fuerza-desplazamiento, para la diagonal 
de acero del contraviento de la estructura de la figura 10, que es de 100/100/10 
mm., y tiene una longitud de 3.80 m., si a = 0.05 y la ductilidad /^^ = 10. Presentar 
el modelo constitutivo para: 

i) El elemento trabajando a compresion. 

ii) El elemento trabajando a traccion. 

En la tabia 3 se indican las propiedades del perfil tubular de 100/100/10 

mm. 



TabIa 3 Propiedades geometricas de un perfil tubular de 1 00/1 00/1 0 mm 



Propiedad geometrica 


Variable 


Valor 


Radio de giro 


r 


3.6968 cm. 


Area 


A 


36 cm2 


Momento de Inercia 


I 


492 cm4 


Momento Estatico 


z 


122 cm3 



• Solucion 

i) Elemento a compresion 
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3.80 T 
= 102.79 ^ = 12 * 703.069 = 8436.8 



r 0.036968 

Pc[,=A* = 0.0036 * 8436.8 = 30.3725 T. 
Pc= As*fy = 0.0036 * 25300 = 91.08 T. 
Py = mm{Pci,Pc) = 30.3725 T. 
EA = 20430000 * 0.0036 = 73548 T 
30 3725 

Acy= = 0.00041296 m 

'^^ 73548 

Ac£/= 10Acy= 0.0041296 m 

Py = Py + aEA (Acy " Acy ) = 30.37 + 0.05 * 73548 * (0.0041296 - 0.00041296) 
= 44.04 T 

ii) Elemento a traccion 

Py = 91.08 T. 

91.08 

^TY= ^^r-.r. = 0.0012 m. 
73548 

Atu= 10 A7-y= 10 * 0.0012 = 0.012 m 

Py = Py + aEA (Aj.y - Aj.y ) = 91.08 + 0.05 * 73548 * (0.012 - 0.0012) = 130.79 T 

En la tabia 4 se resumen los valores encontrados y en la figura 15 se 
presentan los modelos constitutivos, a la izquierda para compresion y a la derecha 
para traccion. 

TabIa 4 Punto notables del diagrama bilineal fuerza-desplazamiento, para 

diagonal de acero. 



Punto Notable 


Compresion 


Traccion 


Fuerza de fluencia (T.) 


30.37 


91.08 


Desplazamiento de fluencia 

(m.) 


0.00041 


0.0012 


Fuerza ultima (T.) 


44.04 


130.79 


Desplazamiento ultimo (m.) 


0.00412 


0.012 
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0 0 5 1 1 5 2 2 5 3 3 5 4 4 5 

Desplazamiento (m) ^^q-^ 




"0 0 002 0 004 0 006 0 008 0 01 0 012 

Desplazamiento (m) 



Figura 16 Modelos constitutivos fuerza desplazamiento para diagonal que 
trabaja a compresion o traccion. 

ELEMENTOS DE ACERO 



En la tabia 5 se presenta la propuesta del ASCE 41 para hallar el 
diagrama momento-rotacion a partir del punto de fluencia; la parte superior 
corresponde a las vigas y la parte inferior a las columnas. En el sistema de 
computacion CEINCI-LAB se programo las secciones indicadas en la figura 16, el 
significado de las variables, en forma general es el siguiente: b^tp es el ancho y 
espesor del patin; h es la altura del perfil; t^^, es el ancho del alma en los perfiles 1 
y 2 0 es el ancho del perfil tubular. 

TabIa 5 Parametros y criterios de aceptacion para procedimientos no lineales 

en vigas y columnas de acero. 



Vigas - Flexion 


A 


b 


c 


bf ^ 52 h ^ 418 

a. — < -= y — < -= 






0.6 


bf ^ 65 h ^ 640 

b. — >-=o — > 


40y 


69 y 


0.2 
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c. Otros 


Interpolacion lineal, se usara el resultado menor. 


Columnas - Flexion 


A 


b 


c 


Para P/PcL < 0.2 








bf ^ 52 h ^ 300 

a. < -j= y — < -j= 


99y 


IWy 


0.6 


, Of ^ DD tl ^ 4oU 

b. — > -= o — > -= 




eOy 


0.2 




Interpolacion lineal, se usara el resultado menor. 


Para 0 2 < P/P^, < 0 ^ 








6f ^ 52 /l ^ 260 

a. < y — < -j= 




[ 5 P \ 


0.6 


, bf ^ 65 h ^ 400 

b. — > -= o — > -= 


Wy 


l.Sdy 


0.2 


c. Otros 


Interpolacion lineal, se usara el resultado menor. 



* Columnas con P/Pa > 0.5 deberan ser controladas por fuerza axial. 




secci6n tipo 1 



bi 



tw 



bs 



tfi tfs 

secci6n tipo 2 



bs 




SECCION TIPO 3 



Figura 17 Secciones de acero consideradas en el sistema CEINCI-LAB 



Para ingresar a la tabia 5 se deben controlar algunas relaciones 
geometricas como — no hay problema con las unidades pero en las relaciones 

donde interviene fy^ se debe colocar la fluencia del acero en ksi. Cuando se 
trabaja con un acero A 36, el valor de fye = 36 ksi. 

La parte inferior de la tabIa 5, corresponde a columnas (flexo compresion) 
y primero se debe controlar la relacion — , donde P es la carga axial que gravita 

PCL 

en la columna y Pr, es la carga critica de pandeo. Cuando la relacion — > 0.5 

PcL 

debe calcularse como elementos sujetos solo a fuerza axial. 

Una vez que se ha clasificado de acuerdo al valor de — , se pasa a 

PcL 

controlar los parametros relacionados con la geometria de la seccion, en forma 
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similar a lo que se hizo con las vigas y nuevamente es probable que se tenga que 
interpolar. 



• Ejempio 5 

El perfil de acero A 36, indicado en la figura 18 es W 30x99 y tiene un 
momento de inercia / = 1660.76 cm*, un momento estatico z = 5122.8 cm^; tiene 
una longitud L = 7.315 m. Se pide: 

i) Encontrar el diagrama momento rotacion, empleando la 
recomendacion del ASCE 41 . Considerar a = 0.03 (relacion de rigidez 
post fluencia con respecto a la elastica). 

ii) Hallar el diagrama momento curvatura para una longitud plastica Lp = 
0.15 m. 



J bf=26.77 cm ^. 



4 tfs=1.7018cm 



h-75.438 +iw=l.3208cm 

4 tfi=1.7018cm 
Figura 18 Perfil W 30x99 de Ejempio 5. 

Solucion 

Relacion momento rotacion 



bf _ 26.77 
2tf~ 2 * 1.7918 

h 75.438 



= 7.6390 



52 



52 
V36 



1.3208 



= 57.1154 



V fye 
300 300 



= 8.6667 



fye V36 



= 50.00 



Al calcular los parametros que se hallan en la tabia 5 se puede determinar que 
para el caso a) cumple con la primera relacion pero no con la segunda. 



52 



2tf 



'f 4fye 

Asi que, probaremos con el caso b) 



pero 



300 



Jye 



65 65 

— = ^= = 10.8333 
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460 460 



fye V36 



= 76.6667 



Para el caso b) no cumple con ninguna relacion por lo que se tendra que 
realizar una interpolacion lineal para encontrar las relaciones. 



2tf 



< 



65 



lye 



h 460 

y —< -j^ 



Tabia 6 


nterpolacion lineal para los parametros A,b y c 


A 


b 


c 






119y 


0.6 


50.00 


7.6659 dy 


9.6659 9y 


0.4933 


57.1154 


^9y 


69y 


0.2 


76.6667 



7.315 * 0.0051128 * 25300 

0v = = 0.0046 rad 

y 6 * 20'430000 * 0.00166076 

Ou = 7.6659 By = 0.0356 rad 
My = z*fy = 0.0051128 * 25300 = 129.9538 T m 
My = My + 0.03 EI (du - By) 
My = 129.9538 + 0.03 * 20'430000 * 0.00166076 (0.0356 - 0.0046) = 161.508 T m 
M^ = cMy = 0.4933 * 129.9538 = 64.1062 T m 
Br = 9.6659 By = 0.0445 rad 
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ii) 




0 015 



0 02 0 025 0 03 0 035 

Rotacion (rad) 



0 0 005 0 01 

Figura 19 Diagrama momento rotacion de perfil W 30x99 
Relacion momento curvatura 



My 

0y = = 



129.9538 rad 

= 0.0038 

EI 20'430000 * 0.00166076 m 



e„ 0.0356 rad 

0(/ = 0y + / = 0.0038 + — — = 0.2411 

Lp 0.15 m 



0 04 



C.045 




0 01 0 016 

Curvatura (rad/m) 



0 021 
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Figura 20 Diagrama momento curvatura bilineal de ejempio 5. 



• Ejempio 6 

En la figura 21 se presenta un perfil W 14x193 de acero A 36 con el cual 
se ha construido una columna que tiene una altura de 3.81 m., sobre la que gravita 
una carga axial de P = 100 T. El momento de inercia / = 99896 cm*; el momento 
estatico z = 5817.4 cm^; un radio de giro r = 0.1651 m. Se pide: 



i) Encontrar la relacion momento rotacion, para un factor a = 0.03 
(relacion entre rigidez post fluencia con respecto a elastica). 

ii) Hallar la relacion momento curvatura, de un modelo bilineal, para una 
longitud plastica Lp = 0.15 m. 



.. bf=39.878 cm .. 

* if 



h=39.37 



= tfs=3.6576 cm 



tw=2.2606 cm 



II tfi=3.6576cm 

Figura 21 Perfil W 14x193 de ejempio 6. 

• Solucion 



Relacion momento rotacion 

L 3.810 T 
- = = 23.0769 ^ F^r = 12 * 703.069 = 8444.3 

r 0.1651 ^ 

PcL=A* F„ = 0.036645 * 8444.3 = 309.4414 T. 

P 100.00 

= 0.3222 



Pci 309.4414 



P 

0.2 < — < 0.5 



bf 39.878 52 52 

= ^ ^^r-r,. = 5.4514 -— = -— = 8.6667 

2tf 2 *3.6576 R- V36 



ye 

h 39.37 260 260 

— = ——= 17.4157 -= = -= = 43.333 

tw 2.2606 IJ- V36 
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Al calcular los parametros que se hallan en la tabia 5 se puede determinar que 
cumple para el caso a), para lo cual se tiene las siguientes expresiones para los 
parametros a.byc 

bf 52 h 3600 
- — < pero — < 

4fye -Jfye 

/ 5 P \ I 100 \ 

« = ^H^ - 3 d 'y = - 309:^^1^) 'y = 'y 

/ 5 P \ I 100 \ 

^ = - 3 d = " ^^^) = 

c = 0.60 

LzL/ P\ 3.81 * 0.0058174 * 25300 / 100 \ 

0v = ^ 1 = 1 = 0.0031 rad 

y 6EI \ PcJ 6 * 20'430000 * 0.00099896 V 309.4414/ 

Bu = 5.0753 * 0.0031 = 0.0157 rad 

Br = 7.8437 * 0.0031 = 0.0243 rad 

MY = z*fy = 0.0058174 * 25300 = 147.1802 T m 

My = My + 0.03 EI (By - By) 

Mfj = 147.1802 + 0.03 * 20'430000 * 0.00099896 (0.0157 - 0.0243) 
= 154.8947 
M,f = cMY = 0.6 * 147.1802 = 88.3081 T m 



Bjt = 9.6659 By = 0.0445 rad 
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ii) Relacion momento curvatura 



My 147.1802 rad 

= — = = 0.0072 

^ EI 20'430000 * 0.00099896 m 

e,, 0.0157 rad 

0(7 = 0y + / = 0.0072 + — — = 0.1119 

Lp 0.15 m 
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7. DIAGONALES RIGIDIZADORAS DE ACERO 

En la tabia 7 se presentan los valores de a,b,c, con los cuales se halla el 
diagrama momento-rotacion de diagonales rigidizadoras, de acero, que se 
encuentran trabajando a compresion o a traccion. En este caso el criterio de 
control viene dado por la longitud efectiva de pandeo kl para el caso de las 
diagonales que trabajan a compresion. 

Donde K es el factor de longitud efectiva y r es el radio de giro del perfil. El 
procedimiento se debe verificar si el elemento es compacto o no para lo cual 
existen las primeras expresiones. 



TabIa 7 Parametros y criterios de aceptacion para procedimientos no lineales 
de diagonales rigidizadoras de acero. 



Diagonal Rigidizadora en 
Compresion 


A 


b 


c 


a. Para Secciones Esbeltas 








Kl/r > A.2^E/fy 








1. Secciones W,I,2L,2C en 
piano 


O.SAc 


lOAc 


0.3 


2. 2L,2C fuera de piano 


O.SAc 


9Ac 


0.3 


3. HSS, Tubos redondos y 


10.5 


79 


0.3 


rectangulares 


b. Para Secciones compactas 
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Kl/r < 2.1^E/fy 








1. Secciones W,I,2L,2C en 
piano 




8Ac 


0.5 


2. 2L,2C fuera de piano 


lAc 


7Ac 


0.5 


3. HSS, Tubos redondos y 
rectangulares 


lAc 


7Ac 


0.5 


Diagonal Rigidizadora en 
Tension 




UAr 


0.8 


Vigas, Columnas en tension 


5At 


7At 


1.0 



8. ESTRUCTURA DE ACERO CON DISIPADORES DE ENERGIA 

En la figura 24 se presenta una estructura de acero, de seis pisos, que 
tiene luces iguales de 7.315 m., y altura de entrepisos, tambien iguales de 3.81 m., 
con disipadores de energia sobre contravientos en forma de "V" invertida. En el 
ejempio 5 se analizo una viga W 30x99 y en el ejempio 6, una columna W 14x193. 
La estructura se analizo con disipadores TADAS, que se muestran a la izquierda 
de la figura 25 o ADAS que estan a la derecha, de la mencionada figura 23. 
(Aguiar etal. 2015, 2). 
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Figura 25 Descripcion de los disipadores de energia. 



En la figura 26 se muestra a la izquierda la numeracion de nudos y 
elementos. A la derecha en cambio se indican los grades de libertad, en este case 
se tiene 84 grades de libertad de los cuales los 12 primeros son coordenadas 
principales. 
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En la figura 27 se presentan los diagramas momento curvatura de algunos 
elementos estructurales, hallados con un valor a - 0.03 y un valor = 10. 
Cuando se encuentra la respuesta en el tiempo, con un analisis no lineal, en cada 
instante de tiempo se obtiene el diagrama momento-rotacion o el diagrama 
momento-curvatura y para el caso de las columnas, estos depende de la carga 
axial, de tal manera que los presentados son referenciales y tienen como objetivo 
que el lector pueda reproducirlos utilizando el sistema de computacion, CEINCI- 
LAB, que se indica en el siguiente apartado pero unicamente para hallar los 
diagramas indicados utilizando las tablas proporcionadas por el ASCE 41 . 




Cur«turai,rad/'m , Cup«tura(rad/m) 



Figura 28 Diagrama momento curvatura de elemento disipador, en la cara 

superior. 

En la figura 28 se presenta los diagramas momento curvatura del 
elemento disipador. Notese que los disipadores de los tres primeros pisos, son 
iguales y el espesor es de 3 cm., y el numero de placas es n = 10. Los disipadores 
del cuarto piso son de 4 mm, existe una reduccion considerable y el del ultimo piso 
es de 2 mm. Por este motive es que disminuye considerablemente la capacidad a 
flexion de los disipadores 64 a 66. 

• Ejempio 7 

Encontrar el diagrama momento-curvatura de un disipador TADAS en el 
borde superior, si el ancho es de 15 cm., el espesor (altura) es de 3 cm., el numero 
de placas es n = 10 (Elemento 61). El material es acero A 36; considerar a = 0.03 
y Me = 10. 

• Solucion 



b = 0.15 m 



h = 0.03 * 10 = 0.3 m. 
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0.15*0.3^ , , 0.3 

I = = 3.375 * 10-4 ^4 — ^ 0.15 m. 

12 2 
/ 3.375 * 10-* 

z = - = — — = 0.0022 

c 0.15 

My = z/y = 0.0022 * 25300 = 56.925 Tm 
EI = 20430000 * 3.375 * 10"* = 6895.1 Tm. 

My 56.925 rad 

0Y = —= = 0.0083 

^ EI 6895.1 m 

rad 

0c = 10 * 0y = 0.083 

m 

My = 56.925 + 0.03 * 6895.1 * (0.083 - 0.0083) = 72.377 Tm 
Los resultados, corresponden al grafico del elemento 61 de la figura 27 



1 0. COMENTARIOS Y CONCLUSIONES 



El objetivo principal de este articulo, es ensenar al lector, en forma 
sencilla, la forma de obtener los diagramas momento-rotacion, utilizando los 
criterios y tablas del ASCE 41, y como se pasa al diagrama momento-curvatura, 
por intermedio de la longitud plastica. Para lograr este objetivo se presentaron 
algunas deducciones como el calculo de la rotacion de fluencia, la longitud 
plastica, entre otras, y se realizaron algunos ejemplos en vigas y columnas de 
hormigon armado y acero; asi como en una diagonal de acero se obtuvo la 
relacion fuerza-desplazamiento, tanto a traccion como a compresion. 

Por otra parte, se ilustro la forma como se obtiene el diagrama memento 
curvatura para un disipador de energia ADAS o TADAS, toda vez que las 
estructuras con estos dispositivos de control son adecuadas para zonas sismicas. 

Finalmente se han presentado programas del sistema de computacion 
CEINCI-LAB, que permiten calcular los diagramas momento-rotacion; momento- 
curvatura; fuerza-desplazamiento, siguiendo los lineamientos del ASCE 41. 
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